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LONG-TERM  GOALS 

My  long-term  goal  is  to  develop  numerical  codes  which  are  applicable  to  the  full 
range  of  low-frequency  surfzone  motions,  to  test  and  refine  these  codes  using 
results  from  held  measurement  programs,  and  to  apply  these  codes  to  the  pre¬ 
diction  and  characterization  of  fluid  motions  and  resulting  sediment  transport 
and  evolution  of  beach  morphology. 

OBJECTIVES 

The  scientific  objective  of  this  study  is  to  gain  an  understanding  of  various 
features  of  the  low  frequency  nearshore  current  climate,  including: 

1.  The  long  term  evolution  and  dynamics  of  shear  waves. 

2.  The  effect  of  interaction  and  competition  between  several  unstable  wavenum¬ 
ber  modes  of  finite  amplitude  shear  waves. 

3.  The  dynamics  of  edge  waves  in  an  environment  where  the  underlying  long¬ 
shore  current  structure  is  strongly  modified  by  the  presence  of  shear  waves. 

4.  The  effect  of  organized  forcing  of  the  shear  waves  by  the  incident  waves. 

The  motions  being  studied  can  provide  a  sizeable  portion  of  the  total  energy 
content  in  inner  surfzone  motions,  and  thus  have  a  (presently  not  understood) 
impact  on  sediment  transport  and  the  evolution  of  beach  morphology.  We  seek 
to  eventually  address  these  issues  using  the  information  about  the  hydrodynamic 
environment  obtained  from  the  present  study. 

APPROACH 

We  have  developed  a  numerical  code  with  the  objective  of  studying  the  motions 
and  combinations  of  motions  described  above.  This  code  has  been  used  to  study 
the  long-time  evolution  of  shear  waves  over  realistic  beach  topographies  with  and 
without  longshore  bars.  We  have  investigated  the  relative  importance  of  bottom 
friction  and  turbulent  and  dispersive  lateral  mixing  mechanisms  in  controlling 
the  growth  of  shear  waves,  and  we  have,  in  turn,  investigated  the  relative  impor¬ 
tance  of  small-scale  lateral  mixing  and  shear  wave-induced  Reynolds  stresses  in 
determining  the  distribution  of  wave-induced  longshore  current.  The  code  will 
be  further  utilized  to  investigate  the  effects  on  shear  wave  evolution  caused  by 
spatial  and  temporal  inhomogeneity  of  both  driving  forces  and  bottom  topog¬ 
raphy.  Finally,  the  model  is  being  extended  to  allow  for  the  computation  of  an 
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evolving  bottom. 

WORK  COMPLETED 

A  model  of  the  nonlinear  shallow  water  equations  with  additional  bottom  friction 
and  forcing  terms  has  been  developed.  The  model  is  robust  and  the  numerical 
code  has  been  refined  in  a  number  of  ways  to  pursue  numerical  experiments 
simulating  various  types  of  low  frequency  motions  in  the  surf  zone  as  well  as  in 
the  swash  zone.  The  model  has  been  extensively  tested  and  has  been  applied  to 
the  generation  of  subharmonic  edge  waves  and  their  growth  to  finite  amplitude 
(Ozkan-Haller  and  Kirby,  1997a). 

The  model  was  used  to  analyze  the  dynamics  of  shear  waves  in  the  absence  of 
other  low  frequency  motions.  For  this  purpose  we  first  simulated  shear  instabili¬ 
ties  of  the  longshore  current  on  a  plane  beach  (Ozkan  and  Kirby,  1995).  Studies 
were  also  carried  out  examining  if  shear  instabilities  of  the  longshore  current 
account  for  the  low  frequency  energy  observed  at  the  Superduck  experiment 
(Ozkan-Haller  and  Kirby,  1996).  The  low  frequency  climate  at  Superduck  was 
examined  by  attempting  to  identify  the  combination  of  bottom  friction  and  mix¬ 
ing  coefficients  which  led  to  a  flow  regime  that  best  reproduced  the  observations 
(Ozkan-Haller  and  Kirby,  1997b). 

RESULTS 

Studies  of  shear  instabilities  on  a  plane  beach  have  produced  results  which  are  in 
agreement  with  independent  results  by  Allen,  Newberger  and  Holman  (J.  Fluid 
Mech,  1996)  who  utilized  a  rigid  lid  model.  It  was  shown  that  the  long  term 
evolution  of  instabilities  in  the  longshore  current  is  strongly  dominated  by  sub¬ 
harmonic  transitions.  These  transitions  were  analyzed  in  detail  (Ozkan-Haller, 
1997).  We  found  that  the  transitions  occur  in  the  form  of  vortex  collisions  where 
a  weaker  vortex  structure  catches  up  and  collides  with  a  slower  more  energetic 
vortex  in  front  of  it.  During  the  collision  of  the  vortices,  there  appears  to  be 
an  exchange  of  identities  of  trailing  and  leading  events,  much  as  in  the  elastic 
collision  of  shallow  water  solitons.  After  the  collision,  most  of  the  energy  is 
transferred  to  the  initially  weaker  trailing  vortex  and  a  subsequent  reduction 
in  the  number  of  large  vortices  is  evident.  The  resulting  flow  structures  are 
longshore  progressive  and  exhibit  strong  offshore  directed  velocities.  These  re¬ 
sults  are  possibly  suggestive  of  a  mechanism  for  the  formation  of  migrating  rip 
currents.  Since  it  is  suspected  that  rip  current  motions  may  be  intensified  by 
the  distortion  to  the  incident  wave  field  caused  by  wave-current  interaction  ef¬ 
fects,  we  are  exploring  the  inclusion  of  this  distortion  in  the  forcing  terms  in  the 
numerical  model. 

Several  cases  using  measured  bottom  bathymetry  and  wave  heights  from  the 
Superduck  field  study  have  also  been  studied.  Bottom  friction  as  well  as  mo¬ 
mentum  mixing  due  to  turbulence  and  depth  variations  in  the  currents  are  taken 
into  account.  A  friction  coefficient  Cf  and  a  mixing  coefficient  M  control  the 
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nature  of  the  simulated  motions.  We  deduce  a  realistic  range  of  values  for  Cf 
and  M  from  measurements.  Decreasing  Cf  within  that  range  results  in  more 
energetic  shear  instabilities  with  higher  propagation  speeds.  At  present,  we  use 
the  predicted  propagation  speed  of  shear  waves  as  the  criterion  for  determining 
a  “best”  friction  coefficient  (see  Figure  1),  and  then  use  the  predicted  maxi¬ 
mum  longshore  current  velocity  as  a  verification  of  this  procedure.  The  shear 
instabilities  induce  significant  momentum  mixing  in  the  surf  zone  and  alter  the 
mean  longshore  current  profile.  The  current  profile  before  the  instabilities  are 
initiated  along  with  the  final  mean  current  profile  are  shown  in  Figure  2  for 
one  case  at  Superduck.  Also  shown  are  the  measured  mean  currents.  It  is  ob¬ 
served  that  the  final  current  profile  is  not  altered  by  the  amount  of  mixing  due 
to  processes  other  than  the  instabilities.  If  more  mixing  due  to  turbulence  is 
assumed  (higher  M),  the  instabilities  are  less  energetic  and  the  mixing  due  to 
the  instabilities  decreases  proportionally  to  produce  the  same  current  profile.  In 
all  cases  mixing  due  to  the  instabilities  is  seen  to  dominate  over  other  mixing 
mechanisms.  However,  the  reason  for  the  apparent  balance  between  the  various 
mixing  mechanisms  and  the  resulting  insensitive  longshore  current  profile  is  not 
understood  and  deserves  further  investigation. 

Vortex  interactions  of  the  type  observed  for  the  simpler  case  of  a  sloping  bottom 
are  also  present  for  the  cases  involving  the  barred  Superduck  bathymetry.  Snap¬ 
shots  of  the  vorticity  field  (Figure  3)  show  the  complicated  nature  of  the  motions. 
Vortex  pairs  propagate  in  the  longshore  direction  y,  interact,  occasionally  merge 
and  are  shed  offshore. 
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Figure  1:  Frequency-cyclic,  longshore  wavenumber  spectra  S(f,k)  (m3/s) 
longshore  velocities  from  measurements  on  October  16.  The  equation  for 
dispersion  line  (dashed)  is  noted. 
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IMPACT  /APPLICATIONS 

This  study  has  shed  light  on  the  processes  that  are  important  in  the  low  fre¬ 
quency  range  of  the  energy  spectrum,  such  as  interactions  between  low  frequency 
waves  and  response  of  the  low  frequency  environment  to  external  forcing.  As 
shear  waves  have  an  important  effect  on  the  velocity  structure  around  the  breaker 
line  they  are  also  a  plausible  explanation  for  offshore  mixing.  As  a  result  of  this 
study  the  magnitude  of  this  effect  can  be  assessed  and  the  general  understand¬ 
ing  of  nearshore  circulation  and  mixing  patterns  can  be  advanced.  This  study 
can  also  serve  as  a  benchmark  for  other  studies  that  are  not  restricted  to  low 
frequency  motions,  but  where  the  low  frequency  motions  are  embedded  in  higher 
frequency  oscillations,  making  the  processes  difficult  to  identify. 

TRANSITIONS 

The  model  code  developed  here  is  being  used  in  an  ONR-funded  study  of  the 
inversion  of  SAR  images  of  the  water  surface  to  determine  bottom  bathymetry 
(with  R.  A.  Dalrymple,  Delaware),  in  conjunction  with  a  much  more  numerically 
intensive  Boussinesq  wave  model.  The  model  and  results  also  form  the  basis  of 
further  proposed  work  on  edge  wave  and  shear  wave  interactions  (Putrevu,  Kirby 
and  Oltman-Shay,  NSF)  and  on  the  characterization  of  nearshore  motion  using 
coherent  structure  analysis  (Kirby,  ONR). 
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Cf  =  0.0035  and  M  =  0  (solid),  M  =  0.25  (dashed),  M  =  0.5  (dash-dotted)  and 
sled  data  (o).  Also  shown  is  the  mean  current  in  the  absence  of  shear  instabilities 
V  for  Cf  =  0.0035  and  M  =  0.5  (dotted). 
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Figure  3:  Snapshots  of  contour  plots  of  vorticity  q  (1/s)  at  5  minute  intervals 
for  cf  =  0.0035  and  M= 0.25. 
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